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ABSTRACT 

Using the science verification data of the Dark Energy Survey (DES) for a new sample of 106 X-Ray se¬ 
lected clusters and groups, we study the stellar mass growth of Bright Central Galaxies (BCGs) since redshift 
1.2. Compared with the expectation in a semi-analytical model applied to the Millennium Simulation, the 
observed BCGs become under-massive/under-luminous with decreasing redshift. We incorporate the uncer¬ 
tainties associated with cluster mass, redshift, and BCG stellar mass measurements into analysis of a redshift- 
dependent BCG-cluster mass relation, m* cx ( i sx'To^Mq )°’^'*^°'°^(1 +^)~o d±o. 34^ compare the observed 
relation to the model prediction. We estimate the average growth rate since z = 1.0 for BCGs hosted by clusters 
of M 2 oo,z = 10*®'®Mq, at z = 1.0: m, 5 CG appears to have grown by 0.13 ±0.11 dex, in tension at ^ 2.5cr signifi¬ 
cance level with the 0.40 dex growth rate expected from the semi-analytic model. We show that the buildup of 
extended intra-cluster light after z = 1.0 may alleviate this tension in BCG growth rates. 

Subject headings: galaxies: evolution - galaxies: clusters: general - galaxies: groups: general 
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1. INTRODUCTION 

Bright central galaxies (BCGs) are the luminous elliptical 
galaxies residing at the centers of galaxy clusters or groups. 
Once commonly referred to as the brightest cluster galaxies, 
the name bright central galaxy better reflects their special 
nature as the central galaxy of a massive halo. BCGs are 
surrounded by a subsidiary population of satellite galaxies. 
Their centrality and large size sets them apart from the gen¬ 
eral galaxy population. 

Early attention about BCGs started with studies about cD- 
type galaxies, since many BCGs are enveloped by extended 
stellar halos (Matthews et al. 1964). Statements that this pop¬ 
ulation is not consistent with being statistically drawn from 
the global galaxy luminosity function led Tremaine &. Rich- 
stone (1977) to argue that BCGs require a special formation 
process. Analytical and early numerical estimates of their 
growth through dynamical friction and resultant cannibalism 
of cluster galaxies was soon identified as a viable process (Os- 
triker & Tremaine 1975; White 1976; Hausman & Ostriker 
1978; Richstone & Malumuth 1983). Early N-body simula¬ 
tions of merging pairs and groups of galaxies led Dubinski 
(1998) to perform the first N-body study of BCG formation in 
a massive halo formed within a cold dark matter (CDM) cos¬ 
mology. In that study, growth through early merging of a few 
massive galaxies dominated over late-time accretion of many 
smaller systems. 

The modern context of BCG assembly through hierarchical 
growth within an evolving spatial network of dark matter ha¬ 
los is now well established, but detailed understanding of vari¬ 
ous competing astrophysical processes remains elusive. Mod¬ 
els in which BCGs accrete their stellar mass through “dry” 
merging with red and old galaxies produce scaling behavior 
and light profiles in fairly good agreement with observations 
(e.g., Ruszkowski & Springel 2009; Laporte et al. 2013). 

Pure N-body models of dry merging ignore intra-cluster gas 
processes such as cooling and subsequent accretion and star 
formation of baryons onto the BCG. Semi-analytical models 
And that such cooling needs to be mitigated by heating, and 
AGN feedback in a so-called “radio mode” is proposed as the 
solution (Croton et al. 2006; De Lucia & Blaizot 2007). Sim¬ 
ulations with explicit hydrodynamic treatment of the baryons 
are struggling to develop sub-grid models that capture the full 
complexity of the baryon behavior {e.g., Martizzi et al. 2012; 
Ragone-Eigueroa et al. 2013; Martizzi et al. 2014; Pike et al. 
2014). 

While BCG in situ star formation is almost certainly sup¬ 
pressed by the quenching effect of AGN (active galactic nu¬ 
clei) feedback (Eabian 1994, 2012), observational studies 
have found that residual star formation of ^ 10- lOOMQyr"* 
exists in many nearby BCGs (Eraser-McKelvie et al. 2014; 
Liu et al. 2012; Groenewald & Loubser 2014). A most 
puzzling study has observed a BCG starburst of 740 ± 
16OM0yr“' in the z = 0.596 Phoenix cluster (McDonald et al. 
2012). Such a large star formation rate would contribute sig¬ 
nificantly to BCG stellar mass even if it lasted for just 1 Gyr. 

Recent arguments based on local cooling-to-dynamical 
timescales tie together this rich phenomenology in a self- 
regulated precipitation model (Voit et al. 2015, and references 
therein). Idealized hydrodynamic simulations (Li & Bryan 
2014a,b; Meece et al. 2015) support an episodic picture in 
which gas below a cooling threshold (roughly Cooi/Gyn < 10) 
feeds black hole accretion and local star formation, with 
AGN feedback serving as the rectifier that shuts down cool¬ 


ing and allows the cycle to refresh. With HST observations of 
BCGs in the CLASH sample, Donahue et al. (2015) offer ev¬ 
idence that ultraviolet morphologies and star-formation rates 
of BCGs in CLASH clusters display features remarkably sim¬ 
ilar to those anticipated by these simulations. 

The semi-analytical expectations of BCG growth have been 
called into question by a number of observations that re¬ 
port significantly slower build-up of stellar mass over time 
(Whiley et al. 2008; Collins et al. 2009; Lidman et al. 2012; 
Lin et al. 2013b,b; Oliva-Altamirano et al. 2014; Inagaki et al. 
2015). This tension highlights limitations in our current un¬ 
derstanding of BCG formation and motivates the work in this 
paper. 

The production of intra-cluster light (ICL) is another im¬ 
portant process affecting BCG formation over time. The ICL 
contains stars that got dispersed into intracluster space from 
BCGs or BCG mergers (see; Contini et al. 2014). Simulation 
and observational studies show that ICL can make up 5-50% 
of the total cluster/group stellar content (Zibetti et al. 2005; 
Krick et al. 2006; Krick &. Bernstein 2007; Gonzalez et al. 
2007; Toledo et al. 2011; Guennou et al. 2012; Burke et al. 
2012; Montes & Trujillo 2014; Giallongo et al. 2014; Pre- 
sotto et al. 2014; Burke et al. 2015). Details of how the ICL 
is formed and how its properties might vary from cluster to 
cluster remain unsettled (Monaco et al. 2006; Conroy et al. 
2007; Puchwein et al. 2010; Rudick et al. 2011; Cui et al. 
2014; Contini et al. 2014; DeMaio et al. 2015). 

To advance our understanding about the above processes 
and BCG formation in general, it is important that we con¬ 
tinue to refine our measurements of BCG growth. Most 
up-to-date observations are yielding perplexing or even con¬ 
tradictory results on this subject, perhaps because of in¬ 
comparability in their processing BCG observables (MandeT 
baum et al. 2005; Lauer et al. 2007; Bernard! et al. 2007). Lor 
instance, a few studies based on high redshift (z > 1.0) X-ray 
selected clusters (Stott et al. 2010, 2011; Collins et al. 2009; 
Whiley et al. 2008) finds no sign of BCG stellar mass growth, 
while others based on clusters at low and high redshifts do ob¬ 
serve the change (Brough et al. 2002; Lin et al. 2013b; Lidman 
et al. 2012, which included some of the samples from the fore- 
mentioned X-ray studies). On the other hand, deriving BCG 
luminosity and hence BCG stellar mass from imaging data is 
not straightforward, and inconsistent measurements may have 
affected many previous findings about BCG formation. Li- 
nally, BCG mass is known to be correlated with cluster mass, 
which needs to accounted for when studying the change of 
BCG mass over time (see for example, Lidman et al. 2012). 
Advances in our understanding of the nature of the growth of 
BCGs require a careful accounting of all of the ingredients, 
including their measurement uncertainties. 

In this paper, we investigate BCG stellar mass growth us¬ 
ing DES Science Verification (DES SV) data, and a new 
sample of 106 X-ray selected clusters and groups from the 
DES XMM Cluster Survey (XCS), an XMM-Newton archival 
discovery project. Through using this X-ray selected sam¬ 
ple, selection effect on studying BCG’s optical properties are 
greatly alleviated; X-ray selected clusters display a wider 
variety of optical properties compared to optically selected 
clusters (Harrison et al. 2012). The cluster and group sam¬ 
ple spans a redshift range of [0,1.2], and a mass range of 
[3 X IO'^Mq, 2 X IO'^Mq]. While most previous studies on 
this redshift range or cluster mass range are combining dif¬ 
ferent samples or different imaging data sets, we study a sin- 
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Fig. 1.— The mass and redshift distribution of the DES-XCS sample (black 
stars) compai'ed to that of Lidman et al. (2012, red circles). 

gle cluster sample with the deep optical data from DES. In 
this paper, we also pay particular attention to possible biases 
affecting BCG photometry, and have carefully evaluated the 
uncertainties associated with cluster mass, redshift, BCG lu¬ 
minosity and BCG stellar mass measurements. We provide 
details of our uncertainty, bias and covariance estimations in 
Appendices A to D. 

The rest of this paper is organized in the following order. In 
Section 2, we present our data sets and derive cluster masses, 
BCG luminosities and BCG stellar masses. We perform a 
matching exercise of BCG redshift evolution to the Millen¬ 
nium Simulation expectations in Section 3, then fit both sim¬ 
ulated and observed BCG populations to a simple low-order 
model in Section 4. We compare this model to previous esti¬ 
mates of BCG growth rate in Section 5. We summarize our 
results in Section 6. Appendices A to D describe the un¬ 
certainties, biases and covariances of relevant measurements. 
Throughout this paper, we assume to be 0.3, Ha to be 0.7, 
and the Hubble parameter h to be 0.7. 

2. DATA 

This paper is based on an X-ray selected cluster and group 
sample from the DES-XCS project. BCG photometry is de¬ 
rived from DES Science Verification data. 

The rest of this section introduces the DES-XCS sample 
and the DES S V data, and also summarizes our procedures of 
deriving cluster masses, selecting BCGs and measuring BCG 
properties. Appendices A to C should be considered as exten¬ 
sions of this section. 

2.1. DES Science Verification Data 

The Dark Energy Survey is a ground-based optical survey 
that uses the wide-held DECam camera (Elaugher et al. 2015) 
mounted on the 4m Blanco telescope to image 5,000 deg^ of 
the southern hemisphere sky (Sanchez 2010). The paper is 
based on 200 deg^ DES Science Verihcation (SV) data. This 
data set was taken during the 2012B observing season be¬ 
fore the main survey (Diehl et al. 2014) began. A large frac¬ 
tion of the SV data have full DES imaging depth (Lin et al. 
2013a) and are processed with the official DES data process¬ 
ing pipeline (Mohr et al. 2012). A more detailed review can 
be found in Sanchez et al. (2014). 

2.2. The DES-XCS Cluster and Group Sample 


The XMM Cluster Survey serendipitously searches for 
galaxy cluster (and group) candidates in the XMM-Newton 
archive (Lloyd-Davies et al. 2011; Mehrtens et al. 2012; Viana 
et al. 2013). The cluster candidates are then verihed with op¬ 
tical/infrared imaging data, which confirm the existence of 
red sequence galaxies. Photometric redshifts of the confirmed 
clusters are also subsequently derived with the red sequence 
locus. Using DES SV data. Miller et al. (in prep.; referred 
to as M15 in the rest of the paper) have identified ^ 170 X- 
ray selected clusters and groups from XCS. M15 also mea¬ 
sures their photometric redshifts and verify the measurements 
against archival spectroscopic redshifts In this paper, we 
use a sub-sample from M15 that consists of 106 clusters and 
groups with mass above 3.0 x IO'^Mq. These clusters and 
groups are all referred to as "clusters" in the rest of the paper. 
In Eigure 1, we show their mass and redshift distribution. Eor 
comparison, we also show the mass and redshift distribution 
of the cluster sample used in a similar study (Lidman et al. 
2012). Our sample covers a lower mass range, and appears to 
be more evenly distributed in the redshift-mass space. 

The cluster mass {M 200 , the mass inside a 3D aperture 
within which the averaged matter density is 200 times the 
critical density) is either derived with X-ray temperature or 
X-ray luminosity, using a lensing calibrated M-T relation 
(Kettula et al. 2013). Because XCS is a serendipitous survey, 
not all the clusters have high quality X-ray temperature mea¬ 
surements. Eor these clusters, we derive their masses from 
X-ray luminosity. Eurther details about this procedure, and 
about the mass uncertainties can be found in Appendix A. 

We note that a handful of the clusters do not seem to 
have significant galaxy over-density associated with them. 
It is possible that our sample contains spurious clusters 
which originates from foreground/background X-ray contam¬ 
inations. We have re-analyzed our analysis after removing 8 
clusters that are not associated with significant galaxy over¬ 
density. The results are consistent with those presented in this 
paper within 0.5 a. Given that these 8 clusters are in the low 
mass range (generally below 10'^ Mq), removing them may 
introduce an artificial mass selection effect. We therefore do 
not attempt to do so in this paper. We also note that other 
factors, including BCG photometry measurement and cluster 
mass scaling relations at the low mass end (see discussion in 
Section 4.2 and 5.2), have bigger effect on our results than the 
possible spurious clusters in the sample. 

2.3. BCG Selection 

The BCGs are selected through visually examining the DES 
optical images, the X-ray emission contours, and the galaxy 
color-magnitude diagram. In this procedure, we aim to se¬ 
lect a bright, extended, elliptical galaxy close to the X-ray 
emission center, which also roughly lies on the cluster red se¬ 
quence. If there exist several red, equally bright and extended 
ellipticals close to the X-ray center, we select the nearest one. 
We did not notice a proper BCG candidate with a blue color. 

We check our visual BCG selection against the central 
galaxy choice of a preliminary version of the DES SV 
RedMaPPer cluster catalog (see the algorithm in Rykoff et al. 
2014). Out of the 106 XCS clusters and groups, 64 are 
matched to RedMaPPer clusters and the majority (61) iden¬ 
tify the same BCG. In the cases where we disagree with the 
BCG, we choose the brighter, more extended galaxy closest 

http: / /ned. ipac. caltech. edu 
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Separation Between BCGs and X-Ray centers [Mpc] 


Fig. 2.— Distances between the BCGs and the X-ray emission centers 
for our cluster sample. Flalf of the BCGs are separated less than 0.07 Mpc 
(transverse comoving distance with negligible uncertainties from redshifts 
measurements) from the X-ray centers, and the lai'ge separations (> 0.4 Mpc) 
happen in clusters that may not be relaxed or appear to have spurious fore¬ 
ground/background emissions. 

to the X-ray center while RedMaPPer selects a galaxy further 
away. The other 42 non-matches are caused by the different 
data coverage, redshift limit, and mass selection of the two 
catalogs: the RedMaPPer catalog employs only a subset of 
the SV data to achieve relatively uniform depth for selecting 
rich clusters below redshift 0.9. 

In Figure 2, we show the distance distribution between the 
selected BCGs and the X-ray emission centers. Half of the 
BCGs are separated by less than 0.07 Mpc (comparable to Lin 
& Mohr 2004) from the X-ray emission centers, regardless of 
the redshifts of the clusters. 

2.4. BCG Photometry, Luminosity and Stellar Mass 

Measuring BCG photometry is among the most controver¬ 
sial topics in BCG studies. In Appendix B, we discuss compli¬ 
cations and possible biases associated with Petrosian magni¬ 
tude, Kron magnitude, prohle htting magnitude and aperture 
magnitude with extended details. We use magnitude mea¬ 
sured with circular apertures of 15 kpc, 32 kpc, 50 kpc and 
60 kpc radii. The main results are derived with the 32 kpc 
radius apertures, considering the BCG half light radius mea¬ 
surements in Stott et al. (2011). Detailed rationalization about 
this choice and description about our measurement procedure 
can also be found in Appendix B. 

We correct for galactic extinction using the stellar locus re¬ 
gression method (High et al. 2009; Kelly et al. 2014; Rykoff 
et al. in prep.), and compute BCG luminosities and stellar 
masses using the stellar population modeling technique. We 
employ a Chabrier (2003) Initial Mass Function (IMF) and the 
Conroy et al. (2009); Conroy & Gunn (2010) simple stellar 
population (SSP) models to construct stellar population tem¬ 
plates, and select templates according to BCG DBS g,r,i,z 
photometry. We use the best-ht model to compute the K- 
correction factor and the mass-to-light ratio. We evaluate un¬ 
certainties associated with BCG apparent magnitude, redshift, 
and BCG mass-to-light ratio. Further details about these pro¬ 
cedures can be found in Appendix C. 

3. SIMULATION MATCHING ANALYSIS 

We hrst inspect the redshift evolution of BCG luminos¬ 
ity and stellar mass through matching our data with a semi- 
analytical simulation. We compare BCG luminosities and 
stellar masses to the corresponding values in the simulation. 



Cluster Redshift 

Fig. 3.— Mass and redshift distribution of the cluster sample used for this 
paper and distribution of the simulation clusters drawn from DL07. The black 
data points show the masses and mass uncertainties of the XCS clusters. The 
red data points show the median masses and the 0.158 and 0.842 percentiles 
of the simulation clusters. For clarification, we show the mass distribution of 
the resampled DL07 clusters at the redshift of the corresponding XCS cluster 
with a small offset. 

with diagrams analogous to those presented in many previous 
studies (Collins et al. 2009; Lu et al. 2014; Lidman et al. 2012; 
Tonini et al. 2012; Lin et al. 2013b; Oliva-Altamirano et al. 
2014) that overlay redshift evolutions of the observed and 
simulated BCG properties. The simulation involved in this 
comparison is the De Lucia & Blaizot (2007) semi-analytical 
(SAM) simulation (referred as DL07 hereafter) based on the 
Millennium project (Springel et al. 2005; Guo et al. 2013). 

3.1. Simulation Sample Selection 

Since BCG luminosity and stellar mass are known to be 
correlated with cluster mass, the comparison between obser¬ 
vation and simulation need to be made between clusters of 
similar masses. For each BCG in our sample, we compare it 
to a simulation subsample of 100 BCGs hosted by clusters of 
similar masses and redshifts. The simulation data are selected 
with the following procedure. 

1. Identify simulation clusters with redshifts closest to that 
of the XCS cluster. Ideally, we would have identi- 
hed a cluster sub-sample with their redshift distribution 
matching the redshift uncertainty of the XCS cluster, 
but this is not possible since simulations are stored at 
discrete redshifts. 

2. Select from the redshift sub-sample of 100 clusters with 
their posterior mass distribution (log-normal) matching 
the mass uncertainty of the XCS cluster. Note that we 
are not using the cluster mass function as a prior. Ap¬ 
plication of this prior leads to sampling clusters ^0.1 
dex less massive, but leave the conclusions unchanged. 

Note that in the above procedure, we are not considering 
additional cluster properties beyond M 200 and redshift. There 
is emerging evidence that X-ray selected clusters may be bi¬ 
ased in terms of cluster concentration distribution (Rasia et al. 
2013), but it is un-clear how the bias would affect BCG for¬ 
mation study. We also do not consider the Eddington bias 
associated with Lx- The M 200 of the lowest LxIPx systems 
are derived with Tx- Future studies yielding higher precision 
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Fig. 4.— Compai'ison between the observed and simulated BCG properties. 
We show the median and 0.158,0.842 percentiles of the differences. The dots 
are data points outside the range of the uncertainty whiskers, and the red di¬ 
amonds are moving medians. The blue bands show the fitted linear model 
with dependence on look-back time, encompassing Icr uncertainties. We no¬ 
tice that the obseiwed BCGs are becoming under-massive/under-luminous at 
decreasing redshift. Note that we have not included stellar mass or luminos¬ 
ity measurement uncertainties in the linear fit. We also use a 95% confidence 
interval clipping technique to remove outliers. 

on BCG growth may wish to take these selection effects into 
consideration. 

In Figure 3, we show the redshift and the mass distribution 
of the XCS clusters together with the re-sampled DL07 sim¬ 
ulation clusters. The above procedure produces a simulation 
sub-sample that well resembles the probability distribution of 
the XCS sample. 

3.2. Redshift Evolution of the Observed BCGs 

We directly compute the relative luminosity and stellar 
mass difference between the observed and simulated BCGs, 
as shown in Figure 4. 

We notice that the differences between the observed and 
simulated BCGs change with redshift. The effect suggests 

We are comparing the observer frame DES z band luminosity to the 
observer frame SDSS z band luminosity in DL07. The response curves of 
the DES z band and the SDSS z band are similar enough, that the magnitude 
measurements for one object in the two systems shall be close within 0.05 
mag. We have tested this statement through cross matching galaxies in the 
SDSS stripe 82 database and the DES Year 1 coadd database. Although it 
is possible to transform between DES z band magnitudes and SDSS z band 
magnitudes, we avoid doing so because the transformation inevitably makes 
assumption about BCG SEDs. 


that the observed BCGs do not grow as rapidly as in DL07 - a 
different redshift evolution history in the observation. We fit 
the differences with a linear dependence on lookback time: if 
the redshift evolution of the observed BCGs is consistent with 
that in DL07, the slope of the linear fit shall be 0. This null 
hypothesis is not favored. 

In Figure 4, we show the linear fitting result with blue bands 
which encompass the Icr uncertainties. The luminosity red¬ 
shift evolution in the observation is different from the sim¬ 
ulation with a 2.5 a significance (0.028 ± 0.011). The sig¬ 
nificance from the stellar mass comparison is lower at 1.3 a 
(0.015 ± 0.012), but BCG stellar mass is less certain (recall 
that it requires a choice for the mass-to-light ratio) and there¬ 
fore the result is noisier. 

The redshift evolution difference shows that the observed 
BCGs become increasingly under-massive/under-luminous at 
decreasing redshift compared to DL07 (compare the result 
to Lidman et al. 2012; Lin et al. 2013b; Oliva-Altamirano 
et al. 2014). At the lowest redshift bin (z ^ 0.1) in Figure 4, 
the observed BCGs appear to be 0.1 to 0.2 dex under- 
massive/under-luminous as a result of a different redshift evo¬ 
lution history. 

Arguably, the above statement relies on a fitting function 
connecting the difference between the observed and simu¬ 
lated BCG properties to redshift. The significance level of this 
statement depends on the exact form of the fitting function. In 
Section 4 and 5, we present stronger evidence on this state¬ 
ment, through modeling the BCG redshift evolution for both 
observational data and simulation data, testing the model and 
eventually showing the model constraints being different in 
the observation and in the simulation. 

In addition, we are not considering BCG luminosity and 
stellar mass uncertainties in this section (they are not included 
in the linear fitting procedure). We also address this in Section 
4 and 5. 

3.3. High Redshift BCGs 

At z > 0.9, we notice that two of the four BCGs in our 
sample appear to be massive/luminous outliers by ^ 0.5 dex, 
which matches previous findings about massive BCGs at z > 
1.0. In Collins et al. (2009), five 1.2 < z < 1.5 BCGs are iden¬ 
tified to be 0.5 ^ 0.7 dex more massive than the DL07 simu¬ 
lation BCGs, and in Liu et al. (2013), a massive z = 1.096 cD 
type galaxy is discovered in a 5 arcmin^ Hubble Deep Field. 
However, after considering cluster mass uncertainty, and the 
BCG luminosity and stellar mass uncertainties, we can only 
detect the over-massive/over-luminous BCG effect with ^ 1 
a significance. 

4. BCG-CLUSTER MASS RELATION 

To further investigate the growth of BCGs, we turn to mod¬ 
eling a redshift-dependent, stellar-to-halo mass relation. We 
refer to this relation as the BCG-Cluster mass relation in this 
paper. Later, in Section 5, we use this relation to model the 
BCG growth rate from z = 1.0 to z = 0. 

4.1. Modeling the BCG-Cluster Mass Relation 

We model the BCG-Cluster Mass Relation as redshift de¬ 
pendent with the following equation, 

■^00 

logm, = logmo -k alog( ~ ) -I- ^log( 1 -I- z). (1) 

X dex = 10' 
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Fig. 5.— Posterior distribution of the parameters, logmo, a, ^ and e, in the BCG-cluster mass relation (Equation 1), based on the BCG stellar mass derived in 
32 kpc aperture. The histogram in each column shows the marginalized distribution of the coiTesponding parameter. Other panels show the correlation between 
two parameters noted on the x and y axes. 


TABLE 1 

Constraints on logmo, o. and a of the BCG-cluster mass relation (Equation 1) and p erom the outlier pruning procedure 



Prior 

15 kpc 

32 kpc 

50 kpc 

60 kpc 

32 kpc (logM 2 oo > 13.85) 

32 kpc (logM 2 oo < 14.70) 

DL()7 

logmo 

[10, 13] 

11.37 ±0.08 

11.52±0.08 

11.60 ±0.09 

11.61 ±0.09 

11.58 ±0.08 

11.49 ±0.08 

11.698 ±0.004 

OL 

[-0.5,0.8] 

0.20 ±0.08 

0.24 ±0.08 

0.30 ±0.08 

0.32 ±0.09 

0.37±0.10 

0.19±0.11 

0.452 ± 0.004 

0 

[-2, 2] 

-0.15 ±0.31 

-0.19 ±0.34 

-0.24 ±0.39 

-0.19 ±0.40 

-0.62 ±0.34 

-0.06 ±0.37 

-0.912 ±0.026 

e 

[0.001,1] 

0.172±0.015 

0.180±0.018 

0.192±0.019 

0.198 ±0.020 

0.169±0.019 

0.186 ±0.020 

0.1628 ±0.0012 

p 

[0.5, 1.0] 

0.970 ±0.017 

0.970 ±0.018 

0.971 ±0.017 

0.970±0.018 

0.965 ±0.019 

0.967 ±0.017 

N.A. 


This equation adopts a power law dependence on cluster 
mass (Oliva-Altamirano et al. 2014; Kravtsov et al. 2014; 
Brough et al. 2008; Moster et al. 2010, 2013) as well as a 
power law dependence on redshift. We choose Mpiv to be 
1.5 X 10'"^Mq, about the median mass of the XCS clusters. 
We also assume that there exists an intrinsic scatter, e, be¬ 
tween the observed BCG stellar mass and this relation, as 
\ogm^: ohs ~ M(log»z»,e^). Hence, the relation contains four 
free parameters; logmo, a, (5 and e. 

We perform a Markov Chain Monte Carlo (MCMC) analy¬ 
sis to sample from the following posterior likelihood; 

logC = -ilog|C| - ^YTc-*Y+logp(Q). (2) 

In this function, Y is a 106 dimension vector (yi,y 2 ---,yi 06 ), 
with the kth element being the difference between the mod¬ 
eled and the observed BCG stellar masses, as; 

(3) 


The covariance matrix, C, in Equation 2 is the combination 
of the covariance matrices for cluster mass measurements, 
BCG stellar mass measurements, redshift measurements and 
the intrinsic scatter. It has the following form; 


C = Cov(m*)-i-a^Cov(logM2oo)+/3^Cov(log(l-i-z))-i-e^I. 

(4) 

Additionally, we implement an outlier pruning procedure 
as we “fit” (or sampling the posterior distribution in Bayesian 
statistics) for the BCG-cluster mass relation, as described in 
Hogg et al. (2010). To summarize this procedure, we adopt a 
set of binary integers Q = (qi,q 2 , ...jqioe) as flags of outliers, 
qjt = 0 indicates an outlier and is correspondingly modified 
as, 

yk = 10gOT,_i: - logOToutlier, (5) 

where momiier is treated as a 5th free parameter. To penalize 
data pruning, we assume a Bernoulli prior distribution for Q, 


yk — ymodel.k yobs,k- 
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characterized by another free parameter p as, 

p(Q) = l[p‘^^a-py~‘‘^- ( 6 ) 

k 

Eventually, the parameters to be sampled from Equation 2 
are logmo, a, /3, e, Q, p, and logmoutUer- More details about de¬ 
riving the posterior likelihood (Equation 2) as well as choos¬ 
ing the covariance matrix can be found in Appendix D. We 
assume uniform truncated priors for all the free parameters 
except Q, and the final result appears to be insensitive to this 
choice. We perform the fitting procedure for both the ob¬ 
served BCGs from the XCS sample and the simulation BCGs 
sampled from the DL07 simulation (Section 3.1). 

4.2. Constraints on the BCG-Cluster Mass Relation 

In Eigure 5, we plot the posterior distribution of logmo, ot, 
/3, e in Equation 1. We also list their marginalized means and 
standard deviations in Table 1. 

The constraint we derive on a agrees well with the re¬ 
ported values from the literature (Oliva-Altamirano et al. 
2014; Kravtsov et al. 2014; Brough et al. 2008). We also no¬ 
tice that a increases with bigger BCG apertures, indicating 
stronger correlation with cluster mass in the BCG outskirts 
(also see Stott et al. 2012). This effect seems to be justi¬ 
fiable, considering an inside-out growth scenario for BCGs 
(van Dokkum et al. 2010; Patel et al. 2013; Bai et al. 2014). 
Eurther analysis with large apertures is limited by the increas¬ 
ing amount of background noise at BCG outskirt, but a larger 
BCG sample may help quantifying the effect. This effect also 
illustrates the importance of understanding BCG photometry 
measurement when deriving BCG-cluster mass relations. 

Our estimation of logmo, the normalization of Equation 1, 
appears to be lower than the corresponding value in DL07 by 
0.1 - 0.2. As logmo is mainly constrained by low redshift 
BCGs, this result is completely consistent with BCGs being 
under-massive at low redshift as discussed in Section 3.2. 

Our estimation of /3, the index of the redshift component 
in Equation 1, also disagrees with the corresponding value in 
DL07. The constraint on /3 derived from the whole cluster 
sample is different from the simulation value at a significance 
level of 2.3 (T. The constraint from our data is closer to 0, sug¬ 
gesting less change of BCG stellar mass with redshift. Note 
that a further, quantitative conclusion should not be drawn. 
Although /3 is the dominant parameter that describes BCG 
redshift evolution in Equation 1, it is not the only one. The 
mass term in Equation 1 also contains information about BCG 
redshift evolution as cluster M 200 evolves with time. A quanti¬ 
tative analysis of BCG redshift evolution is presented in Sec¬ 
tion 5. 

Our constraint on /3 is highly co-variant with logm* (recall 
the bi-variate normal distribution), but the co-variance shall 
not be interpreted as “degeneracy”; a reasonable m* sampled 
from its marginalized posterior distribution does not make /? 
consistent with the simulation. We also notice that different 
conventions for BCG magnitude measurement can bias the 
constraint on /3. Eor example, using the Kron magnitude from 
the popular SExtractor software (Bertin & Arnouts 1996), 
which tends to under-estimate BCG Kron Radius and there¬ 
fore BCG total magnitude (See discussion in Section B.2. 
This effect happens frequently for our intermediate redshift 
BCGs), shifts /? downward by ^ 1 cr. 

We detect hints that the constraints on a and /3 may depend 
on cluster mass (see Table 1). Eor clusters with logM 2 oo above 



Redshift 



Fig. 6. — To derive BCG stellar mass growth rate from BCG-cluster mass 
relations, we will need to derive cluster mass growth history from simula¬ 
tions. In the top panel, we show the halo mass evolution history of ~ 200 
halos with logM 2 oo = 13.8 at ^ = 1.0. We utilize Equation 7 (we are using the 
BCG-cluster mass relation in DL07) to derive the BCG stellar mass growth 
rate shown by the red dashed line. In the bottom panel, we show the test re¬ 
sult for this method (see Section 5.1 for details). Overall, our approach well 
reproduces the average growth rate within \<j for simulation BCGs. 

13.85, we notice stronger correlation between BCG and clus¬ 
ter masses (larger a, compare it to Chiu et al. 2014; van der 
Burg et al. 2014) and steeper redshift evolution (smaller /3) at 
^ 1.0 (T significance level. However, BCGs in low mass clus¬ 
ters (logM 2 oo < 13.85) are possibly over-massive compared to 
our simulation calibrated BCG-cluster mass relation. Evaluat¬ 
ing the masses of low mass clusters and groups through their 
X-ray observables needs to be handled with care. In this pa¬ 
per, we use lensing calibrated M-T relation of galaxy groups 
and clusters to derive M 200 for most low mass clusters (see 
Eigure 10 and Section A). Arguably, the accuracy of X-ray 
inferred masses of low mass clusters is less well characterized 
than the higher mass end. Thus, in our growth rate determina¬ 
tions we show the difference after excluding the lowest mass 
systems (logM 2 oo < 13.85, about 10% of the sample). Eor con¬ 
sistency, we also examine the effect of excluding the highest 
mass systems (logM 2 oo > 14.7, about 10% of the sample). 

5. BCG STELLAR MASS GROWTH SINCE Z ~ 1.0 

In this section, we compute the BCG stellar mass growth 
rate since redshift 1.0. We derive the growth rate conve¬ 
niently using the redshift-dependent BCG-cluster mass rela¬ 
tion from the previous section. Doing so, we are assuming 
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that a redshift-dependent BCG-cluster mass relation not only 
describes the relation between BCG stellar mass and cluster 
mass at various redshifts, but also describes how BCG stel¬ 
lar mass evolves with time. There is no new measurement 
made with observational data in this section. The redshift- 
dependent BCG-cluster mass relation derived in the previous 
section is the only input from observational data. Our method, 
however, do need new input from simulation data, which is the 
mass evolution history of clusters. 

In this section, we compute the stellar mass growth rate 
for the BCGs hosted by clusters of logM 2 oo = 13.8 at z = 1.0. 
The choice is made as the XCS sample well represents these 
clusters and their low redshift descendants (see Figure 6 for 
the mass evolution history of clusters with logM 2 oo = 13.8 at 
z = 1.0). The method is also applied to clusters of different 
masses, but we do not notice significant change of the conclu¬ 
sions. 


5.1. Method and Test 


We need to know how the cluster mass evolves with red- 
shift in our method. To acquire this information, we select a 
sample of halos with z ^ 1.0, logM 2 oo "^13.8 from the Millen¬ 
nium simulation, and extract their evolution history by iden¬ 
tifying descendants of these halos all the way to z = 0 (using 
the descendantid keyword). We then compute the mean M 200 
evolution of these halos, shown in Figure 6. 

The second step is to use the BCG-cluster mass relation to 
derive the average stellar mass of the BCGs hosted by these 
halos at different redshifts. From Equation 1, the average 
BCG stellar mass relative to some normalization epoch, zo, 
can be expressed as: 


log 


m*,zo 


= alog 


M2Q0,z 

M2m,zo 


+ /31og' 


1 -i-z 

1+zo 


(7) 


We take log-^^ from the above equation as describing the 
average BCG stellar mass growth. The component in 

the equation is the average cluster mass growth extracted from 
the simulation. 

The result of applying Equation 7 to the average M 200 
growth with the simulation data is also shown in Eigure 6. 
We estimate the uncertainties on the BCG stellar mass growth 
rate through sampling the joint constraint on a and /3. We do 
not consider the uncertainties of cluster mass evolution as it is 
marginal and is cosmology dependent. 

We test our method by applying it to the DL07 simulation 
BCGs. We first derive the BCG-cluster mass relation in DL07 
using the procedure from Section 4.1 for the sample drawn 
from Section 3.1. We compare the computed BCG growth 
rate to the values obtained through directly tracking cluster 
descendants. The latter is acquired through recording the cen¬ 
tral galaxy stellar mass of the halo descendants since redshift 
1.0. We consider the result from this second approach as the 
true growth of simulation BCGs. 

In the bottom panel of Eigure 6, we show the BCG growth 
rate derived with Equation 1, and the true growth rate en¬ 
compassed by uncertainty from bootstrapping. Overall, for 
low mass clusters, our approach reproduces the average BCG 
growth rate from z = 1.0 to z =0 within Itr. Bias associ¬ 
ated with this method (like progenitor bias, see: Shankar et al. 
2015), if there is any, appears to be negligible. 


5.2. Growth Rate from z = 1.0 to z = 0 



Fig. 7.— This plot shows our BCG stellar mass growth estimation with the 
full sample and with two mass-limited sub-samples. We also show the mea¬ 
surements in Lidman et al. (2012); Lin et al. {2013b) and the BCG growth 
rate in the DL07 simulation. Our estimation is consistent with previous mea¬ 
surements, but slower than DL07 by 2.5 a. The uncertainty from DL07 is 
extremely small as the simulation is well sampled (see Figure 6). 

We compute the BCG stellar mass growth rate using Equa¬ 
tion 7 and compare it to the simulation value obtained with 
the same method. We discuss the observational result based 
on BCG 32 kpc aperture stellar masses in this paper - the re¬ 
sult derived with other apertures look similar. Erom z = 1.0 to 
z = 0, we estimate the BCG growth rate tobe0.13±0.11 dex, 
comparing to 0.40 ±0.05 dex in simulation (uncertainty esti¬ 
mated for the BCG sample in Section 3.1), as shown in Eig¬ 
ure 7. This result is in agreement with our conclusion from the 
simulation matching analysis (Section 3.2), and also in agree¬ 
ment with previous studies (Lidman et al. 2012; Lin et al. 
2013b). Even after considering all the uncertainties, biases 
and covariances associated with BCG luminosity and stellar 
mass measurements, we still confirm that the observed BCG 
growth is slower than the prediction from DL07 at a signifi¬ 
cance level of ~ 2.5 a. 

Like our constraint on the BCG-cluster mass relation, our 
result here shifts by ^ 1 cr (0.29 ± 0.11 dex) when we exclude 
the lowest mass systems (logM 2 oo < 13.85, about 10% of the 
sample). Note that the shift may be caused by inaccuracy of 
X-ray cluster mass scaling relations at the low mass end (see 
discussion in Section 4.2). Eor consistency, we also show the 
result (0.07 ±0.12 dex) after excluding the highest mass sys¬ 
tems (logM 2 oo > 14.7, about 10% of the sample). Our result is 
also susceptible to improper BCG magnitude measurements. 
Using the Kron magnitude from SExtractor, the result will be 
biased toward more rapid BCG growth by ^ 1 cr. We also 
considered applying our method with stellar-to-halo mass re¬ 
lations from literature, but as many previous studies are based 
on magnitude conventions with various problems for BCGs 
(see discussion in Appendix B), we opt for not using them in 
this paper. 


5.3. Role ofICL Production 

In this paper, we have shown from two different perspec¬ 
tives that the BCG stellar mass growth rate in clusters with 
logM 2 oo = 13.8 at z = 1.0 is slower than the prediction naively 
expected in a hierarchical formation scenario (De Lucia & 
Blaizot 2007). This effect is not that surprising with a second 
thought on the processes that contribute to (or counter-act) 
BCG formation. 

A hierarchical structure formation scenario predicts that 
galaxy mergers add stars to BCGs. The BCG stellar build-up 
can be further augmented by in situ star formation, but a re- 










10 



0.0 0.2 0.4 0.6 0.8 1.0 

Redshift 


Fig. 8. — We show BCG growth rates from toy models adding more rapid 
ICL production or more dramatic BCG star fonnation to the DL07 simula¬ 
tion. Including extra constant ICL production at 20 - 40 M^/yx well repro¬ 
duces our observed growth rate. We also show the predicted BCG growth 
rate from Contini et al. (2014), which has updated the DL07 simulation with 
more realistic ICL production. The BCG growth rate in Contini et al. (2014) 
agrees well with our measurement. 

duction in stellar mass is possible from mergers that eject stars 
into the intracluster space. The competition between these 
mechanisms remains a subject to large modeling uncertain¬ 
ties in simulations. If we assume that BCGs experience the 
rapid build-up events (mostly merging events) as prescribed 
in the DL07 simulation, there must be a mechanism that off¬ 
sets BCG growth to mimic slower evolution we observe in 
this paper. 

In Figure 8 , we experiment with incorporating extra stellar 
mass gain or loss into the DL07 simulation. Stellar mass gain 
tends to steepen BCG growth over time, while stellar mass 
loss tends to slow down BCG build-up and flatten the BCG 
growth curve. In order to explain the observed BCG growth 
rate in our data, the BCGs in DL07 would need to go through 
extra stellar mass loss at 20 - 40 Mq/yp ending up with 2.0- 
3.5 X IO'^Mq at z = 0, which agrees with our data. 

Such stellar content stripping from mergers would produce 
intra-cluster light (ICL). Our result indicates that ICL accu¬ 
mulates at 20-40 M 0 /yr after z = 1.0, totaling (1.5-3) x 
1 O"M 0 ICL in the present epoch. This amount corresponds 
to about 30% - 60 % the total of BCG and ICL stellar masses, 
consistent with the observed ICL fraction in low and medium 
redshift clusters (z < 0.5: Zibetti et al. 2005; Krick et al. 
2006; Krick & Bernstein 2007; Gonzalez et al. 2007; Toledo 
et al. 2011; Montes & Trujillo 2014; Giallongo et al. 2014; 
Presotto et al. 2014). 

In fact, ICL production has already been suggested as an ex¬ 
planation to the seemingly mild evolution of massive galaxies 
(Monaco et al. 2006; Conroy et al. 2007; Burke et al. 2012; 
Behroozi et al. 2013; Oliva-Altamirano et al. 2014). Although 
not completely settled, recent studies indicate that ICL possi¬ 
bly forms late, mostly after z=l.0 (Contini et al. 2014; Conroy 
et al. 2007). Specifically, the Contini et al. (2014) study up¬ 
dates the DL07 simulation with more realistic ICL production 
processes, and predicts slower BCG growth rate (Figure 8 ), 
in excellent agreement with our measurement. Hence, the 
slow BCG stellar mass growth since z = 1.0 observed through¬ 
out this paper is completely justifiable if ICL forms late after 
z=L0. 

Admittedly, the DL07 simulation also includes stellar strip¬ 
ping that would produce ICL. Unfortunately the amount of 
ICL from this simulation is not retrievable, and we are not 
able to analyze if it meets our expectation. The Guo et al. 


(2011) SAM simulation has explicitly included ICL produc¬ 
tion and predicts very similar BCG growth with DL07, but 
much of the ICL is already in place before z = 1.0, which is 
not favored in our interpretation. 

6. SUMMARY AND DISCUSSION 

Using new photometric data from DBS and a new X-ray se¬ 
lected cluster and group sample from the XCS, we investigate 
the redshift evolution of BCG stellar mass since z = 1.2. We 
derive constraints on the BCG-cluster mass relation, and com¬ 
pute the BCG stellar mass growth rate for our sample. From 
two different perspectives, we demonstrate that the BCG stel¬ 
lar mass growth since z = 1.0 is slower than the expectation 
from a semi-analytical simulation implementing a simple hi¬ 
erarchical BCG formation scenario. The discrepancy is de¬ 
tected with a significance level as high as 2.5 a. We find this 
slow growth rate after z = 1.0 to be compatible with the late 
formation of ICL. 

We have carefully considered various uncertainties related 
to studying BCG growth in this work, including the uncertain¬ 
ties of BCG stellar mass measurements, cluster/BCG redshift 
measurements and cluster mass measurements. We explic¬ 
itly consider these uncertainties through likelihood analysis, 
and expect this analysis to help clarify ongoing discussions 
about how statistical and systematic uncertainties affecting 
BCG growth measurements. 

We also adopt a simple but novel method to compute BCG 
stellar mass growth rate. Despite considerable attention paid 
to this topic in the literature, BCG stellar mass growth has 
been studied with various techniques inconsistent with each 
other. Ideally, one would like to evaluate BCG stellar mass 
growth by comparing the BCG masses within the same clus¬ 
ter at high and low redshifts, as we did for method testing in 
Section 5.1. This is not possible with observations. However, 
Lin et al. (2013b) have adopted the idea through constructing 
a cluster sample that resembles the average halo evolution his¬ 
tory. In observational studies, the more common approach is 
to compare the BCG masses of a high redshift cluster sample 
and a low redshift cluster sample, while adjusting the cluster 
mass binning at different redshifts to account for cluster mass 
evolution (Collins et al. 2009; Lidman et al. 2012; Bai et al. 
2014). The results from these observational studies are widely 
compared to De Lucia & Blaizot (2007), which computes the 
BCG stellar mass growth rate through a “fixed space density” 
method, i.e, selecting the 125 most massive clusters at z ^ 1.0 
and z ~ 0 respectively to compare their BCG masses. Com¬ 
pared to these previous studies, our method allows consistent 
comparison to simulation for clusters of specific masses and 
redshifts. Our test in Section 5.1 shows that the approach suf¬ 
fers from only negligible bias for the required precision. 

Finally, the analyses presented in this paper are based on 
DBS SV data, a data set corresponding to only 5% of the nom¬ 
inal DBS footprint. With spectroscopic and X-ray follow-up. 
Miller et al. (in prep) show that the final DBS/XCS sample 
should be about 10 times larger than this data set. Comparing 
the constraints on the BCG-cluster mass relation derived with 
1000 simulation clusters rather than 100 of them, we conclude 
that we expect ^ 3 times improvement in the measurement 
uncertainty of BCG growth. At this level of statistical power, 
it will be critically important to thoroughly understand the un¬ 
certainties associated with various observables. This paper 
presents the first steps toward such an analysis. 
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APPENDIX 

A. GALAXY CLUSTER MASS 
A.l. Cluster Mass From X-ray Temperature 

We use a lensing calibrated M—T scaling relation from Kettula et al. (2013) to derive cluster mass from X-ray temperature 
{Tx, core not excised). In Kettula et al. (2013), weak lensing mass measurements are obtained for 10 galaxy groups in the mass 
range 0.3-6.0 x IO'^/iyqMq. Together with 55 galaxy clusters (many above 2 x from Hoekstra et al. (2011, 2012); 

Mahdavi et al. (2013), Kettula et al. (2013) derive weak lensing calibrated M—T relation across the group to cluster range 
We check the Kettula et al. (2013) M-T scaling relation against a few other studies based on the gas content (Figure 9). The 
Hydrostatic Equilibrium (HSE, Sun et al. 2009; Eckmiller et al. 2011; Vikhlinin et al. 2009) and gas mass fraction (Mantz et al. 
2010) calibrated relations agree with the Kettula et al. (2013) M- T relation at the cluster scale, but have troubles matching to it 

Chandra temperatures in Mahdavi et al. (2013) are adjusted to match 
XMM calibration. 
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Fig. 9.— Comparison of the cluster masses derived from a few M —T scaling relations. We plot the mass derived from the Kettula et al. (2013) relation on the 
a:— axis of all the panels. The Kettula et al. (2013) M-T scaling relation agrees well with other relations at the cluster scale (See Section A.l). For simplification, 
we only include X-ray temperature measurement uncertainty in this figure. 


at the group scale (Sun et al. 2009; Eckmiller et al. 2011). As known from simulations (Nagai et al. 2007; Rasia et al. 2012, see 
Kettula et al. (2013) for discussion), the disagreement is not surprising as HSE masses are biased low at the group scale and the 
gas mass fraction relation is only derived with the most massive clusters. 

We use the Kettula et al. (2013) scaling relation to derive cluster M 500 from X-ray temperature, and then the Hu & Kravtsov 
(2003) relation to derive M 200 from M 500 . We assume the cluster concentration parameter to be 5 in this procedure. Using a 
different concentration parameter in the [3, 5] range only changes M 200 at percent level. We also assume the intrinsic scatter of 
M 500 to be 0.1 dex as typically found in simulation studies (Kravtsov et al. 2006; Henry et al. 2009). 

A.2. Cluster mass from X-ray Luminosity 

We resort to X-ray luminosity (Lx, core not excised) to estimate masses for clusters/groups that do not have high quality 
temperature measurements. We first derive X-ray temperature using a L-T scaling relation, and then derive M 200 using the 
procedure above. 

We use a self-similar, redshift-dependent L-T scaling relation from Hilton et al. (2012), but also experimented with a few 
other self-similar L— T relations (see Eigure 10) from Stott et al. (2012); Maughan et al. (2012); Pratt et al. (2009). The Hilton 
et al. (2012) relation, which is also based on a XCS sample (Mehrtens et al. 2012), provides the best fit to our data. We assume 
0.1 dex intrinsic scatter for the derived temperatures, as it is constrained in Hilton et al. (2012). 

A.3. Mass Uncertainties and the Choice between Temperature and Luminosity Based Masses 

We decide between Lx and Tx based masses through comparing their uncertainties. To estimate the mass uncertainties as¬ 
sociated with each method, we produce 200 "pseudo-measurements" for each cluster, sampling through temperature/luminosity 
measurement uncertainty, the scaling relation uncertainty, and the intrinsic scatter of the relations. We derive mass uncertainties 
for Lx or Tx based masses assuming log-normal distribution for the 200 "pseudo-measurements". If the uncertainty of Tx mass 
is larger than the uncertainty of Lx mass by 0.05 dex (we prefer Tx mass since Lx mass is more susceptible to biases), we use Lx 
mass in lieu of Tx mass. In the end, about half of the cluster masses are derived with Lx, and most of the clusters masses have < 
0.25 dex uncertainty. 

The 200 "pseudo-measurements" are also used to derive the M 200 covariance between the cluster sample. Because we are 
including scaling relation uncertainty, the M 200 covariance matrix is not diagnonal. 
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Fig. 10.— Left: We derive cluster masses from both X-ray luminosities and X-ray temperatures, and decide which one to use through compai'ing their 
uncertainties. Right: X-ray temperature and X-ray luminosity of the XCS clusters plotted against the scaling relations in Hilton et al. (2012); Stott et al. (2012); 
Maughan et al. (2012); Pratt et al. (2009).The grey band shows the redshift-dependent L-T relation in Hilton et al. (2012) between z = 1 and z =0. 


B. BCG PHOTOMETRY 

BCG photometry measurement is based on products from the official DBS Data Management (DESDM) processing pipeline 
(Mohr et al. 2012). In this pipeline, single exposure images are processed, calibrated and later background-subtracted and co¬ 
added to produce coadd images. DESDM also runs an advanced version of SExtractor (Bertin 2011; Bertin & Arnouts 1996) on 
processed single exposure images and coadd images to produce star/galaxy catalogs, which we do not use because of existing 
problems for BCGs. In this paper, we derive BCG photometry from processed single exposure images. The following Sections 
B. 1 to B.4 describe our explorations on measuring BCG flux with different magnitude conventions. We discuss potential problems 
associated with Petrosian magnitude, Kron magnitude, profile fitting magnitude, and aperture magnitude here, but most of the 
problems are already well analyzed in literature (especially, see Graham &. Driver 2005; Haussler et al. 2007; von der Linden 
et al. 2007; Bernard! et al. 2014). Nevertheless, we provide a summary in this Section. 

Our final choice is to use BCG photometry measured with 15 kpc, 32 kpc, 50 kpc, and 60 kpc apertures for this paper. 

B.l. Petrosian Magnitude 

Petrosian magnitude measures the flux enclosed within a scaled aperture known as the "Petrosian radius", which is calculated 
considering background noise level and object light profile (Petrosian 1976; Blanton et al. 2001; Yasuda et al. 2001). It is 
extraordinarily robust under exposure to exposure variations, but not appropriate for extended galaxies. Although Petrosian 
magnitude accounts for most of the flux of a disk-like (Sersic index =1) galaxy, it will only recover 80% of the flux for a 
bulge-like galaxy with a de Vaucouleurs (Sersic index =4) profile (Blanton et al. 2()01). 

Indeed, a series of studies have found that using Petrosian magnitude (see Bernard! et al. 2013; He et al. 2013, for relevant 
discussion), the brightness of Luminous Red Galaxies (LRGs) is under-estimated by about 0.3 mag. Moreover, the missing flux 
problem is sensitive to the profiles of extended galaxies, and worsens quickly with higher Sersic index. Graham & Driver (2005) 
estimate that Petrosian magnitude at its most popular configuration (one that is adopted by SDSS) under-estimates the luminosity 
of Sersic index =10 galaxies by 44.7% (0.643 mag)! Eor this reason, we are not exploiting Petrosian magnitude in this paper. 

B.2. Kron Magnitude 

Kron magnitude is another scaled aperture magnitude, measuring the flux enclosed within a few "Kron radius" (usually 2.5 
Kron radius), and the Kron radius is decided from the light profile (Kron 1980). Kron magnitude is not as robust as Petrosian 
magnitude under exposure to exposure variations, but does appear to be more proper for extended galaxies. 

Like Petrosian magnitude, Kron magnitude recovers most of the flux of a disk-like galaxy, but misses 10% of the flux for a 
bulge-like galaxy (Sersic index ^ 4, Andreon 2002; Graham & Driver 2005). Unlike Petrosian magnitude, the flux missing ratio 
is in-sensitive to the galaxy Sersic index. Graham & Driver (2005) estimate that the missing flux varies only at percent level 
when Sersic index changes from 2 to 10. Indeed, tests with simulated skies (see Andreon 2002, or our test in Eigure II) show 
that Kron magnitude only underestimates the brightness of buIge-Iike galaxies by about 0.2 mag. It also appears to be indifferent 
to the presence of ICE; when we apply the measurement to simulated BCGs enclosed by ICE (we use the model in Giallongo 
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Fig. 11.— We investigate measurement bias associated with Kron magnitude and aperture magnitude using the UFIG sky simulation (see Section B.4 for 
details). Note that this test is done for a general galaxy population rather than BCGs. BCGs below redshift 1.0 generally have Sersic index > 2 and apparent 
magnitude below 22. In the top row, mime is the galaxy’s input total magnitude, but in the bottom row, mime is the galaxy’s input 32 kpc aperture magnitude. 
To summarize this figure, Kron magnitude tend to under-estimate the brightness of bulge-like galaxies and extended galaxies, while aperture magnitude remain 
well-behaved for galaxies of all profiles and sizes. The measurements from both systems do become biased for faint galaxies with apparent magnitude above 
23, but the bias is un-important for this work. For efficiency, we use SExtractor output in this paper, which compare well with our own measurements (see 
Appendix B.4 for a description on the procedure ) for a general galaxy population. 
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Fig. 12.— The populai' SExtractor software tends to under-estimate BCG Kron radius, resulting in significant brightness under-estimation. In panel (b), we 
show the difference in Kron magnitude measurement. Am, when the measurement is made to a different Kron radius (x-axis). 
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Fig. 13.— (a): Measurement dispersion of aperture magnitude, profile fitting magnitude and Kron magnitude, (b): We sample BCG apparent magnitude from 
multiple exposures to evaluate the measurement uncertainty, and our measurement accuracy is limited by the number of exposures we have. 

et al. 2014), the measurement changes only < 0.1 mag. 

As proper as the design of Kron magnitude seems to be, the real problem comes from observationally deriving the Kron 
radius. As pointed out in Graham & Driver (2005), correctly estimating Kron radius requires integration over the light profile to 
a very large radius, usually many times the half light radius for extended galaxies. If the integration is improperly truncated, the 
measured Kron radius will be much smaller, and Kron magnitude turns out to be catastrophically wrong - it may under-estimate 
the flux of an extended galaxy by as much as 50% (Bernstein et al. 2002)! 

We find this to be a frequent problem for BCG measurements from the widely-used SExtractor software (i.e., mag_auto), as 
demonstrated in Figure 12 (a) and (b). The Kron radius from SExtractor is two times smaller than it should be for one of the 
BCGs, and the BCG light intensity at 2.5 SExtractor Kron radius is still high. As a result, SExtractor underestimates the Kron 
flux of this BCG by ~ 0.5 mag. This problem seems purely algorithmic though. Using the galaxy intensity profile to re-calculate 
Kron radius until it converges, we are able to correct this measurement error. Comparing the corrected measurements to the 
magnitude measurements from profile fitting (see Section B.3), we recover the 0.2 mag accuracy of Kron magnitude as discussed 
above. 

For this paper, we have re-done our analysis using Kron magnitude. We re-compute the Kron radius instead of using SExtractor 
output, but the result remains qualitatively similar. 

B.3. Profile Fitting Magnitude 

We have also experimented with BCG profile fitting magnitude from the GALFIT software (Peng et al. 2002, 2010). We fit the 
BCGs w2015PASP..127.1183Zith a model consisting two Sersic profiles, one with Sersic index = 1 (i.e., a disk profile) and one 
with flexible Sersic index as suggested by Bernardi et al. (2014); Meert et al. (2015). We convolve these models to point spread 
functions (PSF) derived with the PSFex software (Bertin 2011), and carefully mask all neighboring objects including blended 
objects identified with the GAIN deblender (Zhang et al. 2015). Overall, the design of this procedure is similar to the Galapagos 
fitting software (Barden et al. 2012). 

For this paper, we only use the profile fitting magnitude for testing purpose (see Section B.2 and Section B.4). We hesitate 
about using it for scientific purpose as we realize that the measurement needs to be extensively tested with sky simulations as 
in Haussler et al. (2007); Bernardi et al. (2014); Meert et al. (2015). Upon evaluating the profile fitting magnitude uncertainties 
(see Figure 12), we do not find it to improve BCG measurement accuracy and therefore do not consider the testing efforts to be 
worthwhile for this paper. We nevertheless have re-done our analysis using this magnitude, but the result remains qualitatively 
similar. 


B.4. Aperture Magnitude Used in this Paper 

In this paper, we measure BCG magnitude with circular apertures of 15 kpc, 32 kpc, 50 kpc, and 60 kpc radii. The main results 
in this paper are derived with the 32 kpc measurements, considering the BCG half light radius measurements in (Stott et al. 2011). 
The 32 kpc aperture choice is also comparable to the popular Kron magnitude aperture (2.5 Kron radius) measurements from the 
SExtractor software. We carefully mask BCG neighbors (including blended objects identified with the GAIN deblender, Zhang 
et al. 2015) and interpolate for the BCG intensity in the masked area. To realistically evaluate BCG magnitude uncertainty, we 
perform the procedure on processed single exposure images, use the median as the measurement, and evaluate the uncertainty 
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TABLE 2 

Parameters oe the stellar population models 


Formation Redshift 

20,10, k 5,4, 3, 2.5,2,1.5,1.0,0.k 0.5 

Metallicity 

0.03, 0.025, 0.02, 0.015, 0.01, 0.008, 0.005, 0.003, 0.002 

R-toldmg time (Ciyrs) 

30, 15, 10, 8, 5, 3, 2, 1, 0.8, 0.5, 0.3, 0.1 

Observed Redshift 

1.50, 1.49, 1.48, ..., 0.03, 0.02, 0.01 


through bootstrapping. We hnd our typical measurement uncertainty to be ~ 0.4 mag, signihcantly larger than the SExtractor 
estimation from co-added images (but not larger when we bootstrap the SExtractor measurements from single exposure images). 
Since we perform the measurements independently on different exposures, our uncertainty is more comprehensive than the 
SExtractor uncertainty from co-added images (also see the magnitude measurement scatter test in Eigure 11). Our measurement 
becomes uncertain when we have few exposures to work with (see Eigure 12), which will be improved as DES assembles more 
exposures in the coming years. 

To evaluate the sky background level around BCGs, we use background check maps generated with the SExtractor software 
from DESDM, conhgured with the "Global evaluation" process. We sample the values in a ring with inner and outer radius of 
~ 13 arcsec and 18 arcsec from the BCG. We have investigated how sky background estimation affects our measurement, as it 
was considered a difficult task for BCGs. We hnd it to have only marginal inhuence. Indeed, even by using a "Global" setting, 
SExtractor still overestimates the background around some extremely bright sources, known as the dark halo problem within DES 
(after background subtraction, the light intensity of a bright object falls slightly below 0 at the outskirt). However, changing the 
background sampling location only marginally shift our hnal measurements. In fact, other details of the measurement procedure, 
like in-complete masking of neighboring sources may cause bigger problems. 

We test for measurement bias associated with aperture magnitude and Kron magnitude with simulated DES images, using the 
UEIG simulation (Berge et al. 2013; Chang et al. 2014). This sky simulation is based on an N-body dark matter simulation popu¬ 
lated with galaxies using the Adding Density Determined GAlaxies to Lightcone Simulations (ADDGALS) algorithm (Wechsler 
in prep.; Busha in prep.; Dietrich et al. 2014, for a review). We hnd that both aperture magnitude and Kron magnitude tend to 
under-estimate the brightness of fainter sources, but the effect is negligible for even the furthest BCGs (z band apparent magnitude 
is about 22). In addition, the bias would only have suppressed the signihcance of our result, as further objects are evaluated to be 
less massive/luminous. We also perform the test with sky simulations based on adding simulated galaxies into real DES co-add 
images, known as the Balrog simulation (Suchyta et al. 2015), and came to the same conclusion. 

C. BCG STELLAR MASS UNCERTAINTY 

We derive BCG luminosity and stellar mass with the stellar population modeling technique, and use a SED htting procedure to 
hnd a stellar population model for each BCG. This procedure begins with using the EZGal package (Mancone & Gonzalez 2012), 
the Chabrier (2003) Initial Mass Eunction (IME) and the Conroy et al. (2009); Conroy & Gunn (2010) simple stellar population 
(SSP) models to produce stellar population templates with various star formation histories and metallicities. We make templates 
with exponentially-decaying star formation histories (the r model) characterized by various e-folding time, metallicity, formation 
redshift and observed redshift. In Table 2, we list the parameter values used for these templates. 

We then use a Chi-Square minimizing technique (see: Mitchell et al. 2013) to decide the stellar population template for each 
BCG. The htting procedure is done with BCG photometry in DES g, r, i, z bands and we ht the BCGs only to templates of their 
observed redshifts. After a best ht is identihed for each BCG, we use the K-correction term from the template to compute BCG 
luminosity, and then the mass-to-light ratio to compute BCG stellar mass. We derive BCG luminosity in DES z band, and BCG 
stellar mass according to z band luminosity. As an alternative, we also use the Blanton & Roweis (2007) K-correction package to 
derive BCG luminosity, but the result remains unchanged. 

C.l. BCG Luminosity Uncertainty 

We estimate BCG luminosity uncertainty combining BCG magnitude and redshift (see M15) uncertainties. To simplify sub¬ 
sequent analyses, we assume the redshift uncertainty to be statistical (systematic uncertainty is about ^ 0.001, comparing to ^ 
0.05 statistical uncertainty, see M15). The redshift uncertainty is taken as 0.001 if archival spectroscopic redshift is available. We 
ignore K-correction uncertainty as it is well decided. 

C.2. BCG Mass-to-Light Ratio Uncertainty 

We estimate BCG stellar mass uncertainty combining BCG luminosity uncertainty and BCG mass-to-light ratio (MLR) uncer¬ 
tainty. In this section, we pay special attention to estimating the MLR uncertainty from modeling star formation histories, which 
is the uncertainty from htting a r-model to BCGs that formed through merging with galaxies of various star formation histories. 
We evaluate the uncertainty by applying the stellar population htting procedure to DL07 BCGs (selected in Section 3.1). We 
compare the derived BCG MLR to their true values in the simulation. 

In the left panel of Eigure 14, we show the difference between the derived and true MLR plotted against redshift. We notice 
a systematic uncertainty of ~ 0.05 dex, likely caused by the mismatch of IME in our procedure and in DL07 (a Chabrier IME 
produces a mass-to-light ratio 0.05 dex higher than that of a Kroupa IME, see Papovich et al. 2011). 

We also notice a statistical scatter with the derived values, ranging from 0.05 dex to 0.1 dex with weak dependence on redshift 
(Eigure 14), but no dependence on cluster mass or BCG stellar mass. We evaluate the uncertainty and covariance for our BCG 
sample taking the corresponding values in simulation. To elaborate, for each BCG in our sample, we assume its MLR to have 
been measured 100 times (each BCG is matched to 100 simulation BCGs in Section 3.1), and the error of each measurement is 
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Fig. 14. — (a) We apply our stellar model fitting procedure to the DL07 simulation BCGs using their simulated SDSS g, r, i, z photometiy. Because the fitting 
procedure and the simulated BCG photometry are based on different SSP models and IMFs, there exists a 0.05 dex systematic offset. The derived mass-to-light 
ratio also has statistical uncertainty of 0.05 dex at z = 0 and 0.1 dex at z =1. (b) We evaluate the stellar mass uncertainty combining redshift, magnitude and 
mass-to-light ratio uncertainties. The stellai* mass uncertainty is dominated by apparent magnitude uncertainty at z < 0.8, and redshift uncertainty at z > 0.8. 


the offset between the derived and true MLR for one simulation BCG. As a result of this set-up, the MLR uncertainty for each 
BCG contains about 0.05 dex systematic uncertainty and 0.05 to 0.1 dex statistical uncertainty depending on its redshift. 

Admittedly, it is more than likely that we are underestimating the BCG MLR uncertainty. In our simulation test, the systematic 
uncertainty originates from using slightly different SSP models and IMFs (Conroy et al. (2009); Conroy & Gunn (2010) SSP 
models and Chabrier (2003) IMF in our procedure VS the Bruzual & Chariot (2003) SSP models and Kroupa (2001) IMF in 
DL07). The statistical uncertainty originates from matching t star formation history and fixed metallicity to DL07 BCGs. We 
have not considered uncertainties associated with SSP models, dust distributions and possible IMF variations 
Estimating the uncertainties from these so-called "known unknowns" is difficult. Conroy et al. (2009) shows that one may at 
best recover the MLR of bright red galaxies with 0.15 dex uncertainty at z = 0, or 0.3 dex at z = 2.0. According to this result, we 
would have under-estimated BCG MLR uncertainty by ~ 0.1 dex. We also experimented with the SSP models from Maraston 
(2005) and Bruzual & Chariot (2003), but the derived MLR differences are lower than 0.1 dex. 

Since the redshift dependence of our estimation is qualitatively similar to that presented in Conroy et al. (2009), it is un¬ 
likely that we are affected with our conclusion about BCG redshift evolution. We therefore do not attempt to include additional 
uncertainties from the "known unknowns". Eventually, the BCG stellar mass uncertainty is dominated by the uncertainty from 
magnitude measurement or redshift (See Eigure 14b), rather than from MLR. 


D. COVARIANCE AND THE LIKELIHOOD FUNCTION 

This appendix provides extra information about the likelihood function presented in Section 4. We assume multivariate normal 
distribution for Y. Combined with a Bernoulli prior distribution for Q, the posterior likelihood becomes: 


£ocp(Q)x|C| '^exp(^-^Y^C 'y). 

Re-write the likelihood at the log scale and ignore the normalization of the probability distribution, we have: 

log C = -ilog|C| - iYTc-‘Y + log/7(Q). 


Because 


y = y,,,„,-y 


model 5 


The covariance matrix used in the posterior likelihood can be further expanded as : 

C = Cov(Y,Y) 

~ Cov(Yjnodel~ ^obs? ^model Yobs) 

“ Cov(Yobs; Yobs) “t" Cov(Y£oodeL Ymodel) 

= Cov(m*,obs, ni,,obs) + a^Cov(logM 2 oo, logM2oo)+/3^Cov(log( 1 -i- z), log( 1 -i- z)) -i- ct^I. 


(Dl) 


(D2) 

(D3) 


(D4) 
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We have simplified this expression as Equation 4 in Section 4. 

If the covariance matrices for cluster redshift, cluster M 200 and BCG stellar mass are all diagonal, our posterior likelihood 
function - Equation 2 - would have the same form as that derived in Hogg et al. (2010). In our case, the covariance matrices for 
M 200 and BCG stellar masses are not diagonal because of systematic uncertainties (See details in Sections A to C). 

E. TABLE - BCGS AND THEIR HOST CLUSTER PROPERTIES 
TABLE 3 


Cluster Name ^200 Cluster Redshift BCG Stellar Mass BCG DES z Luminosity 


I 

XMMXCSJ003545.5-431758.6 

14.17 

± 

0.20 

0.650 

± 

0.002 

11.48 

± 

0.07 

-23.45 

± 

0.05 

2 

XMMXCSI003548.1-432232.8 

14.30 

± 

0.16 

0.630 

± 

0.002 

11.53 

± 

0.06 

-23.78 

± 

0.04 

3 

XMMXCSJ021442.21-043315.3 

14.72 

± 

0.12 

0.142 

± 

0.001 

11.75 

± 

0.06 

-24.13 

± 

0.05 

4 

XMMXCSJ021547.28-045031.4 

13.75 

± 

0.23 

0.955 

± 

0.108 

11.75 

± 

0.15 

-24.16 

± 

0.32 

5 

XMMXCSJ021611.69-041422.8 

14.44 

± 

0.17 

0.154 

± 

0.001 

11.52 

± 

0.07 

-23.77 

± 

0.12 

6 

XMMXCSI021735.7-051315.8 

14.25 

± 

0.21 

0.646 

± 

0.002 

11.34 

± 

0.06 

-23.38 

± 

0.07 

7 

XMMXCSJ021807.88-054557.3 

14.56 

± 

0.20 

0.689 

± 

0.002 

11.29 

± 

0.39 

-23.07 

± 

0.96 

8 

XMMXCSJ021826.0-045942.7 

14.26 

± 

0.20 

0.873 

± 

0.002 

10.20 

± 

0.13 

-21.10 

± 

0.30 

9 

XMMXCSI021833.0-050101.5 

14.19 

± 

0.15 

0.874 

± 

0.002 

11.37 

± 

0.14 

-23.70 

± 

0.30 

10 

XMMXCSI021837.90-054037.0 

13.92 

± 

0.20 

0.317 

± 

0.001 

11.65 

± 

0.08 

-24.10 

± 

0.11 

11 

XMMXCSJ021843.7-053253.3 

13.94 

± 

0.18 

0.381 

± 

0.001 

11.10 

± 

0.09 

-22.76 

± 

0.16 

12 

XMMXCSJ021914.5-045053.2 

13.83 

± 

0.39 

0.222 

± 

0.001 

11.40 

± 

0.08 

-23.13 

± 

0.08 

13 

XMMXCSI021940.78-055043.7 

14.23 

± 

0.21 

0.729 

± 

0.002 

11.30 

± 

0.09 

-23.74 

± 

0.14 

14 

XMMXCS1021944.3-045326.8 

14.24 

± 

0.25 

0.332 

± 

0.001 

10.96 

± 

0.11 

-22.20 

± 

0.23 

15 

XMMXCSI022154.83-054519.0 

14.51 

± 

0.20 

0.258 

± 

0.001 

11.67 

± 

0.09 

-23.63 

± 

0.11 

16 

XMMXCSJ022204.5-043246.3 

14.26 

± 

0.19 

0.315 

± 

0.001 

11.23 

± 

0.08 

-23.08 

± 

0.12 

17 

XMMXCSJ022233.02-045803.5 

14.41 

± 

0.20 

0.839 

± 

0.002 

11.68 

± 

0.07 

-24.14 

± 

0.12 

18 

XMMXCSI022302.5-043619.6 

13.84 

± 

0.21 

1.260 

± 

0.085 

11.79 

± 

0.16 

-24.01 

± 

0.20 

19 

XMMXCSJ022318.3-051209.8 

14.00 

± 

0.19 

0.332 

± 

0.001 

11.15 

± 

0.08 

-23.17 

± 

0.03 

20 

XMMXCSJ022318.37-052707.6 

13.69 

± 

0.14 

0.211 

± 

0.001 

11.50 

± 

0.09 

-23.47 

± 

0.12 

21 

XMMXCSI022347.50-025134.4 

13.86 

± 

0.15 

0.185 

± 

0.007 

11.52 

± 

0.14 

-23.64 

± 

0.31 

22 

XMMXCSI022350.88-053643.9 

14.85 

± 

0.20 

0.500 

± 

0.002 

11.47 

± 

0.12 

-23.74 

± 

0.08 

23 

XMMXCSJ022356.32-030556.8 

14.04 

± 

0.20 

0.300 

± 

0.010 

11.31 

± 

0.10 

-23.18 

± 

0.20 

24 

XMMXCSI022357.6-043519.7 

14.44 

± 

0.19 

0.572 

± 

0.002 

11.24 

± 

0.11 

-23.54 

± 

0.20 

25 

XMMXCSJ022401.1-050542.2 

14.04 

± 

0.19 

0.327 

± 

0.001 

11.43 

± 

0.08 

-23.42 

± 

0.12 

26 

XMMXCSI022433.9-041442.7 

14.04 

± 

0.14 

0.262 

± 

0.001 

11.52 

± 

0.21 

-23.56 

± 

0.49 

27 

XMMXCSI022512.2-062307.9 

14.41 

± 

0.14 

0.204 

± 

0.001 

11.74 

± 

0.08 

-24.18 

± 

0.09 

28 

XMMXCSJ022512.31-053112.3 

13.62 

± 

0.20 

0.294 

± 

0.001 

11.54 

± 

0.10 

-23.45 

± 

0.14 

29 

XMMXCSI022524.7-044043.9 

14.18 

± 

0.13 

0.267 

± 

0.001 

11.60 

± 

0.07 

-23.65 

± 

0.04 

30 

XMMXCSJ022529.8-041432.7 

14.00 

± 

0.12 

0.143 

± 

0.001 

11.27 

± 

0.07 

-23.14 

± 

0.07 

31 

XMMXCSJ022532.42-035502.4 

14.02 

± 

0.22 

0.771 

± 

0.002 

11.52 

± 

0.08 

-23.75 

± 

0.09 

32 

XMMXCSI022549.02-055339.3 

14.42 

± 

0.20 

0.232 

± 

0.001 

11.58 

± 

0.09 

-23.55 

± 

0.16 

33 

XMMXCSJ022559.68-024932.4 

14.04 

± 

0.21 

0.305 

± 

0.017 

11.18 

± 

0.15 

-22.91 

± 

0.33 

34 

XMMXCSJ022632.5-054651.9 

14.40 

± 

0.20 

0.445 

± 

0.026 

11.54 

± 

0.10 

-23.86 

± 

0.21 

35 

XMMXCSI022634.75-040408.0 

14.22 

± 

0.21 

0.346 

± 

0.001 

11.67 

± 

0.07 

-24.13 

± 

0.04 

36 

XMMXCSJ022722.82-032147.3 

14.38 

± 

0.13 

0.360 

± 

0.016 

11.57 

± 

0.14 

-23.81 

± 

0.32 

37 

XMMXCSJ022726.4-043206.8 

14.03 

± 

0.25 

0.309 

± 

0.001 

11.58 

± 

0.07 

-23.83 

± 

0.03 

38 

XMMXCSJ022740.4-045129.9 

14.09 

± 

0.12 

0.295 

± 

0.001 

11.37 

± 

0.06 

-23.66 

± 

0.05 

39 

XMMXCSI022802.85-045101.1 

14.57 

± 

0.22 

0.297 

± 

0.001 

11.53 

± 

0.09 

-23.79 

± 

0.06 

40 

XMMXCSJ022808.41-053553.2 

13.98 

± 

0.26 

0.192 

± 

0.001 

11.50 

± 

0.08 

-23.53 

± 

0.06 

41 

XMMXCSI022827.3-042542.5 

14.37 

± 

0.20 

0.433 

± 

0.001 

10.46 

± 

0.15 

-21.66 

± 

0.36 

42 

XMMXCSJ022829.83-031257.2 

14.14 

± 

0.20 

0.370 

± 

0.000 

11.63 

± 

0.13 

-23.98 

± 

0.30 

43 

XMMXCSI023052.5-045128.7 

14.41 

± 

0.20 

0.590 

± 

0.002 

11.34 

± 

0.10 

-25.00 

± 

0.15 

44 

XMMXCSI033151.23-274936.2 

13.77 

± 

0.22 

0.676 

± 

0.002 

11.19 

± 

0.06 

-23.29 

± 

0.04 

45 

XMMXCS J034004.0-283150.6 

14.02 

± 

0.17 

0.262 

± 

0.001 

11.57 

± 

0.06 

-23.47 

± 

0.05 

46 

XMMXCSJ034025.95-284025.4 

14.46 

± 

0.12 

0.067 

± 

0.001 

12.27 

± 

0.06 

-22.38 

± 

0.06 

47 

XMMXCSJ034026.03-282835.8 

14.44 

± 

0.20 

0.375 

± 

0.007 

11.58 

± 

0.07 

-23.93 

± 

0.06 

48 

XMMXCSJ041646.0-552510.4 

14.14 

± 

0.17 

0.410 

± 

0.008 

11.61 

± 

0.10 

-24.08 

± 

0.18 

49 

XMMXCSI042017.8-503155.0 

14.24 

± 

0.11 

0.470 

± 

0.015 

11.60 

± 

0.23 

-24.03 

± 

0.56 

50 

XMMXCSJ042226.36-514029.7 

14.35 

± 

0.13 

0.575 

± 

0.039 

11.69 

± 

0.34 

-24.29 

± 

0.85 

51 

XMMXCSJ043218.04-610356.5 

14.50 

± 

0.23 

0.435 

± 

0.021 

11.67 

± 

0.09 

-24.18 

± 

0.19 

52 

XMMXCSI043318.93-614013.7 

14.55 

± 

0.21 

0.545 

± 

0.007 

11.11 

± 

0.24 

-22.98 

± 

0.51 

53 

XMMXCS J043706.81-541413.0 

14.33 

± 

0.14 

0.505 

± 

0.006 

11.42 

± 

0.20 

-23.56 

± 

0.47 

54 

XMMXCSJ043708.09-542908.8 

14.57 

± 

0.20 

0.565 

± 

0.041 

11.14 

± 

0.16 

-22.91 

± 

0.28 

55 

XMMXCSI043818.09-541917.5 

15.05 

± 

0.13 

0.475 

± 

0.005 

11.84 

± 

0.24 

-24.63 

± 

0.58 

56 

XMMXCS J043940.82-542412.9 

14.67 

± 

0.20 

0.700 

± 

0.015 

11.82 

± 

0.20 

-24.47 

± 

0.47 

57 

XMMXCS J045421.2-531531.1 

13.91 

± 

0.14 

0.325 

± 

0.028 

11.15 

± 

0.11 

-22.49 

± 

0.23 

58 

XMMXCSI045506.1-532343.2 

14.40 

± 

0.18 

0.435 

± 

0.019 

11.51 

± 

0.14 

-23.72 

± 

0.34 

59 

XMMXCS1051141.31-515420.8 

14.44 

± 

0.13 

0.724 

± 

0.028 

11.57 

± 

0.10 

-24.12 

± 

0.22 

60 

XMMXCSJ051636.81-543113.3 

14.84 

± 

0.12 

0.295 

± 

0.001 

12.15 

± 

0.37 

-24.69 

± 

0.92 

61 

XMMXCSI065829.1-555641.5 

15.12 

± 

0.12 

0.297 

± 

0.001 

11.77 

± 

0.16 

-24.10 

± 

0.33 

62 

XMMXCSJ065860.0-560926.8 

14.23 

± 

0.21 

0.335 

± 

0.015 

11.69 

± 

0.08 

-24.13 

± 

0.12 

63 

XMMXCSJ095737.12-1-023426.0 

14.75 

± 

0.18 

0.373 

± 

0.001 

11.99 

± 

0.08 

-24.57 

± 

0.13 

64 

XMMXCSI095823.49-I-024920.2 

14.25 

± 

0.19 

0.343 

± 

0.001 

11.56 

± 

0.08 

-23.64 

± 

0.13 

65 

XMMXCSJ095846.90-I-021550.8 

14.39 

± 

0.19 

0.771 

± 

0.002 

11.38 

± 

0.13 

-23.62 

± 

0.27 

66 

XMMXCSI095924.73-1-014615.7 

13.86 

± 

0.14 

0.124 

± 

0.001 

11.22 

± 

0.10 

-23.02 

± 

0.20 

67 

XMMXCSJ095931.63-1-022657.2 

14.74 

± 

0.21 

0.356 

± 

0.001 

11.57 

± 

0.13 

-23.70 

± 

0.30 

68 

XMMXCSJ095944.68-I-023619.8 

13.53 

± 

0.17 

0.343 

± 

0.001 

11.62 

± 

0.10 

-24.12 

± 

0.19 

69 

XMMXCSI095947.18-1-025529.1 

14.29 

± 

0.19 

0.126 

± 

0.001 

10.95 

± 

0.10 

-22.16 

± 

0.18 
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84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 


XMMXCSJ095951.46+014051.8 

14.25 

± 

TABLE 3 

0.13 0.373 ± 

0.001 

11.69 

± 

0.09 

-24.08 

± 

0.17 

XMMXCSJ100021.72+022329.3 

13.83 

± 

0.15 

0.221 

± 

0.001 

11.40 

± 

0.10 

-23.49 

± 

0.18 

XMMXCSJ100027.26+022135.9 

13.98 

± 

0.21 

0.221 

± 

0.001 

11.33 

± 

0.15 

-23.17 

± 

0.32 

XMMXCSJ100043.28+014608.0 

14.14 

± 

0.24 

0.346 

± 

0.001 

11.58 

± 

0.13 

-23.74 

± 

0.28 

XMMXCSJ100047.16+015917.0 

14.37 

± 

0.11 

0.438 

± 

0.001 

11.45 

± 

0.14 

-23.62 

± 

0.32 

XMMXCSJ100109.18+013336.0 

13.93 

± 

0.20 

0.435 

± 

0.002 

11.64 

± 

0.13 

-24.19 

± 

0.29 

XMMXCSJ100141.74+022538.0 

14.32 

± 

0.21 

0.120 

± 

0.001 

11.58 

± 

0.09 

-23.44 

± 

0.17 

XMMXCSJ100142.56+014059.4 

14.16 

± 

0.21 

0.218 

± 

0.001 

11.25 

± 

0.12 

-23.24 

± 

0.24 

XMMXCSJ100201.42+021334.2 

14.96 

± 

0.21 

0.838 

± 

0.002 

11.57 

± 

0.10 

-23.95 

± 

0.22 

XMMXCSJ232737.63-541610.0 

14.56 

± 

0.21 

1.005 

± 

0.096 

11.43 

± 

0.14 

-24.08 

± 

0.29 

XMMXCSJ232810.21-555015.8 

14.79 

± 

0.19 

0.890 

± 

0.035 

11.62 

± 

0.09 

-24.07 

± 

0.13 

XMMXCSJ232900.4-533931.7 

13.86 

± 

0.17 

0.255 

± 

0.004 

11.68 

± 

0.08 

-23.95 

± 

0.09 

XMMXCSJ233000.57-543706.4 

14.21 

± 

0.12 

0.176 

± 

0.001 

11.67 

± 

0.07 

-23.99 

± 

0.06 

XMMXCS J233003.40-541415.6 

13.94 

± 

0.23 

0.415 

± 

0.011 

11.42 

± 

0.07 

-23.55 

± 

0.09 

XMMXCSJ233037.38-554340.2 

14.09 

± 

0.25 

0.330 

± 

0.017 

11.58 

± 

0.10 

-23.86 

± 

0.18 

XMMXCSJ233135.2-562754.1 

14.42 

± 

0.18 

0.185 

± 

0.005 

11.48 

± 

0.18 

-23.53 

± 

0.41 

XMMXCSJ233204.99-551242.8 

13.82 

± 

0.18 

0.305 

± 

0.014 

11.39 

± 

0.09 

-23.43 

± 

0.14 

XMMXCSJ233215.97-544205.3 

14.21 

± 

0.19 

0.360 

± 

0.022 

11.74 

± 

0.09 

-24.32 

± 

0.18 

XMMXCSJ233331.89-554753.4 

14.53 

± 

0.21 

1.195 

± 

0.065 

11.78 

± 

0.18 

-23.61 

± 

0.33 

XMMXCSJ233346.00-553826.9 

14.53 

± 

0.19 

0.780 

± 

0.000 

11.60 

± 

0.12 

-23.97 

± 

0.26 

XMMXCSJ233406.45-554710.8 

14.46 

± 

0.20 

0.775 

± 

0.000 

11.49 

± 

0.06 

-23.74 

± 

0.12 

XMMXCSJ233421.4-541054.6 

14.22 

± 

0.18 

0.365 

± 

0.017 

11.57 

± 

0.10 

-23.82 

± 

0.19 

XMMXCSJ233429.10-543618.3 

14.42 

± 

0.21 

0.450 

± 

0.009 

11.64 

± 

0.06 

-24.20 

± 

0.10 

XMMXCSJ233528.43-543501.1 

14.63 

± 

0.21 

0.870 

± 

0.118 

11.60 

± 

0.16 

-24.13 

± 

0.37 

XMMXCSJ233637.1-524408.2 

13.90 

± 

0.21 

0.560 

± 

0.011 

11.65 

± 

0.37 

-24.17 

± 

0.93 

XMMXCSJ233706.8-541911.5 

14.23 

± 

0.26 

0.505 

± 

0.008 

11.57 

± 

0.27 

-24.01 

± 

0.67 

XMMXCSJ233745.46-562758.5 

14.27 

± 

0.19 

0.380 

± 

0.018 

11.36 

± 

0.07 

-23.37 

± 

0.11 

XMMXCSJ233836.3-543740.3 

14.52 

± 

0.18 

0.375 

± 

0.006 

11.67 

± 

0.08 

-24.18 

± 

0.10 

XMMXCSJ234119.2-560400.2 

14.39 

± 

0.19 

0.475 

± 

0.014 

11.52 

± 

0.28 

-23.77 

± 

0.70 

XMMXCSJ234142.5-555747.7 

14.37 

± 

0.16 

0.200 

± 

0.005 

11.80 

± 

0.44 

-24.33 

± 

1.09 

XMMXCS J234145.4-564000.7 

14.23 

± 

0.24 

0.495 

± 

0.009 

11.45 

± 

0.09 

-23.77 

± 

0.18 

XMMXCS J234201.2-553253.8 

14.27 

± 

0.21 

0.555 

± 

0.005 

11.26 

± 

0.26 

-23.26 

± 

0.64 

XMMXCS J234231.5-562106.8 

14.30 

± 

0.15 

0.415 

± 

0.027 

11.53 

± 

0.11 

-23.80 

± 

0.22 

XMMXCSJ234448.2-561728.2 

14.30 

± 

0.17 

0.650 

± 

0.006 

11.61 

± 

0.25 

-23.86 

± 

0.59 

XMMXCSJ234730.8-553320.5 

14.05 

± 

0.20 

0.395 

± 

0.025 

11.72 

± 

0.40 

-24.16 

± 

1.00 

XMMXCSJ235009.5-551957.9 

14.08 

± 

0.15 

0.420 

± 

0.015 

11.60 

± 

0.35 

-24.08 

± 

0.85 

XMMXCSJ235059.5-552206.1 

14.20 

± 

0.16 

0.140 

± 

0.006 

11.35 

± 

0.26 

-23.28 

± 

0.63 



